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Abstract 

Hybrid-ARQ protocols have become common in many packet transmission systems due to their 
incorporation in various standards. Hybrid-ARQ combines the normal automatic repeat request (ARQ) 
method with error correction codes to increase reliability and throughput. In this paper, we look at 
improving upon this performance using feedback information from the destination, in particular, using 
a powerful forward error correction (FEC) code in conjunction with a proposed linear feedback code 
for the Rayleigh block fading channels. The new hybrid-ARQ scheme is initially developed for full 
received packet feedback in a point-to-point link. It is then extended to various different multiple-antenna 
scenarios (MISO/MIMO) with varying amounts of packet feedback information. Simulations illustrate 
gains in throughput. 

Index Terms 

hybrid-ARQ, additive Gaussian noise channels, channel output feedback, MIMO fading channel, 
concatenated coding 

I. Introduction 

The tremendous growth in demand for throughput in wireless networks warrants new design principles for cod- 
ing information at the lower layers. In recent years, packet-based hybrid automatic repeat request (ARQ), which 
integrates forward error correction (FEC) coding with the traditional automatic repeat request protocol, has sparked 
much interest. Any hybrid-ARQ scheme includes the transmission of an acknowledgement (ACK) or a negative- 
acknowledgement (NACK) from the destination to the source. Although not normally viewed this way, the feedback 
of ACK/NACK can be seen as a form of channel output information (COI) indicating to the source the 'quality' of 
the channel output. Exploiting the full potential of COI at the source for hybrid-ARQ schemes, however, has not been 
explored in the literature. In fact most of the discussion about feedback in wireless systems has been limited to the use 
of channel state information ( CSI) at the source [ 1 1— [3 1. 

Research since the 1960s [4|-[6| has long established the utility of using COI at the source to increase reliability 
in additive white Gaussian noise (AWGN) channels. The gain in reliability scales very fast in blocklength and relies 
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on simple, linear coding schemes. These schemes achieve a doubly exponential decay in the probability of error as 
a function of the number of packet retransmissions; this is in stark contrast to open loop systems (without COI) that 
can achieve only singly exponential decay in the probability of error. Most of the the literature for COI assumes an 
information-theoretic perspective for analysis. In this work, however, we take a more signal processing approach to 
COI. In particular, we explore the efficacy of including COI to increase the throughput of a practical hybrid-ARQ 
scheme. 

As is commonly done in hybrid-ARQ literature, we assume that the transmissions will take place over a block 
fading channel where the fading characteristics of the channel are assumed to be constant over the transmission of a 
packet (e.g., fl7)— flTD |). In the context of this channel, there will be two main types of side-information that we consider 
to be available at the source: CSI and COI. CSI is commonly known as a complex-valued quantity that represents the 
spatial alignment of the channel. This can be either outdated, where the source has access to outdated values of CSI 
(i.e., the channel states for previous blocks), or current, where the source has access to the present value of CSI (i.e., 
the channel state for the current block). Outdated CSI is commonly obtained through feedback techniques from the 
destination. Current CSI can be obtained through numerous methods such as exploiting channel reciprocity through 
time division duplexing. In fact, most recent standards [11] and technologies like multi-user MIMO, network MIMO, 
and OFDM incur major penalties in performance without the availability of current CSI [ 12 1-[ 14 |. The most common 
form of COI is simply the past received packet at the destination; we will be using this definition as the destination 
incurs no processing before feeding back the received packet. The main focus of the paper will be the integration of 
COI in to a hybrid-ARQ framework. 

Hybrid-ARQ improves the reliability of the transmission link by jointly encoding or decoding the information 
symbols across multiple received packets. Specifically, there are three ways 1031 in which hybrid-ARQ schemes are 
implemented: 

> Type I: Packets are encoded using a fixed-rate FEC code, and both information and parity symbols are sent to 
the destination. In the event that the destination is not able to decode the packet, it rejects (NACK) the current 
transmission and requests the retransmission of the same packet from the source. Subsequent retransmissions 
from the source are merely a repetition of the first transmission. 

> Type IT. In this case, the destination has a buffer to store previous unsuccessfully transmitted packets. The first 
packet sent consists of the FEC code and each subsequent retransmission consists of only the parity bits (in- 
cremental redundancy) to help the receiver at the destination jointly decode across many retransmissions of the 
same packet. 

• Type III: This method is similar to Type II with one major difference. In Type III, every retransmission is self 
decodable, e.g., Chase combining [16]. Therefore the destination has the flexibility to either combine the current 
retransmission with all the previously received retransmissions or use only the current packet for decoding. 



January 19, 2013 



DRAFT 



3 



The first mention of hybrid- ARQ techniques can be traced back to papers from the 1960s (e.g., ifTTl ). However, most 
attention to this protocol has been given during the late 1990s and early 2000s. Throughput and delay analyses were 
done for the Gaussian collision channel in [7|-[9], |[j"8l . These topics were also investigated for wireless multicast in 
1 1 1 and for block fading channels with modulation constraints in |fl9l . The hybrid- ARQ technique has been looked at 
when using many different types of FEC codes including turbo codes Il20l - ll22l . convolutional codes [8), [231 . LDPC 
codes [9 1, [24 1, and Raptor codes [25 1, [26|. The performance analysis of a hybrid- ARQ scheme with CSI [2|, [3|, 
l|27l - l|29l and without CSI at the source [301 has also been studied in the literature. In addition, different ways to 
utilize the feedback channel have been investigated in [ST), [32). 

In this work, however, we explore the advantages of combining conventional CSI feedback with COI feedback. 
As we noted earlier, the potential of COI feedback in hybrid- ARQ has not been properly explored. The method we 
introduce is a variation of Type III hybrid- ARQ scheme that incorporates the use of COI feedback. In the event that no 
COI feedback is available, our proposed scheme simply reduces to the regular Chase combining in which packets are 
repeated for retransmission and the destination combines the received packets using maximum ratio combining (MRC). 
However when COI feedback is available, we look at implementing a linear feedback code that is a generalization of 
Chase combining to increase the performance of the packet transmission system. A linear feedback code is simply 
a transmission scheme in which the transmit value is a strictly linear function of the message to be sent and the 
feedback side-information [33). We show that such codes provide advantages over merely repeating the last packet, 
while offering simpler analysis and implementation than conventional Type II incremental redundancy codes. 

A relevant concern for the implementation of COI feedback techniques is practicality as it requires (possibly) 
sending large amounts of data from the destination back to the source. However, we hope to address these concerns 
by first studying the ideal scenarios (i.e., perfect COI feedback by feeding back the full received packet) to illustrate 
what is theoretically possible and then extend these results to limited-resource cases (i.e., noisy COI feedback and 
feeding back only parts of the received packet). This allows us to establish a trade-off in performance to allow for 
practical limitations on the system. Furthermore, COI feedback techniques can be especially beneficial when there is 
link asymmetry between the source and the destination. In other words, the situations in which the reverse link can 
support much higher rates than the forward link. 

To accommodate the use of multiple-input single-output (MISO) and multiple-input multiple-output (MIMO) 
systems, we first construct the proposed scheme for the simplest case of single-input single-output (SISO) transmission 
and then extend the scheme to the case with multiple transmit antennas. Specifically, the scheme is adapted for use with 
MISO and MIMO when current CSI is available at the source and either perfect or noisy COI feedback is available. It 
is also adapted for MIMO when perfect COI and only outdated CSI is available at the source. 

The paper is structured as follows. In Section UU a brief high-level description of hybrid- ARQ is given to motivate 
the investigation into using more feedback in a packet retransmission scheme. In Section [TTU the feedback scheme 
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to be integrated into a hybrid-ARQ protocol is introduced for SISO systems. We begin the section by describing 
the encoding process. It is followed by a discussion of decoding; this involves two different cases - systems with 
noiseless COI feedback and systems with noisy COI feedback. In Section llVl the SISO scheme is extended to various 
multiple antenna scenarios. In Section [V] the overall hybrid-ARQ system is discussed in detail where now the COI 
feedback schemes created are integrated as a generalization of Chase combining. Schemes that vary the amount of 
COI feedback being sent to the source are also discussed. In Section lVll throughput simulations are given to illustrate 
the performance of the proposed hybrid-ARQ scheme versus other commonly used hybrid-ARQ schemes such as 
||34l and traditional Chase combining. Note that our comparison is with the incremental redundancy ARQ scheme 
in |f34| (not the standardized system in general). It is not our intention to assume the same conditions present in the 
transport channels as the ones discussed in [34|. 



The vectors (matrices) are represented by lower (upper) boldface letters while scalars are represented by lower 
italicized letters. The operators ( ) T , ( )*,tr(-), and ||-|| denote the transpose, conjugate transpose, trace, and Euclidean 
norm of a matrix/vector respectively. The expectation of a random variable or matrix/vector is denoted by E[-}. The 
boldface letter I represents the identity matrix. 



Consider using the SISO hybrid-ARQ transmission system in Fig. Q] where there is one antenna available at the 
source and the destination. The goal of the transmission scheme is to successfully send the binary information packet, 
m 6 GF(2) into , to the destination over a maximum of N packet retransmissions. GF(2) denotes the Galois field 
with just two elements {0, 1}, and £j n f denotes the total number of information bits. Transmission is accomplished 
by first encoding the information packet using a rate iinfo/icodcd (where L; n f < L co d d) FEC code, producing a 
binary codeword of length L C odcd referred to as c € G-F(2) Lcodt!d . The codeword is then modulated using a source 
constellation Q[N] (e.g., QAM, QPSK, etc.) to create a length L packet of modulation symbols called G C lxi . 
Note that the source constellation 9 [N] is a function of the maximum number of transmissions N. If N is large for 
a fixed L (i.e., L/N is small), we may decide to choose a denser constellation to achieve higher throughput. This is 
then processed by a packet encoder that encapsulates most of the hybrid-ARQ process. At this stage, the modulated 
symbols are further encoded to generate the transmitted signal x[k] 6 C lxL . It is worthwhile to contrast that 9 is the 
packet of desired information symbols and the entries of x[fe] are the actual signals sent at each channel use to convey 
that information to the destination. Note that some quantities have a retransmission index, k, which refers to time on 
the packet level (i.e., for each k a length L signal, x[k], is transmitted). Furthermore, the transmit vector is constrained 
by the power constraint at the source given by 



Notation: 



II. System Model 





(l) 
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Fig. 1: The hybrid- ARQ transmission system 



where, as mentioned, N is the maximum number of retransmissions allowed and p > is the average power per 
channel use. 

At the destination, the k th retransmission received signal, y[fc] 6 C lxL , is obtained. Using this setup, y[k] can be 
written as 

y[k] = h[k]x[k] + z[fc], 1 < k < N, (2) 

where z[k] E C lxL is additive noise whose entries are i.i.d. complex Gaussian such that z[k] ~ CJ\f(0, 1), and h[k] g 
C is a zero-mean complex Gaussian random variable with unit variance. Therefore, we assume that the retransmission 
takes place across a Rayleigh block fading channel with blocklength L. Note that x, y, and z have been defined as row 
vectors; this is to aid the later extension of the scheme to a MISO/MIMO setting. After retransmission k, the received 
packet is combined using all previously received packets to create an estimate of the original modulated information 
packet, 0[k]. The combining stage, in Chase combining for example, combines all the received realizations for a given 
symbol using MRC. Improving upon the encoding and combining steps using COI forms the main thrust of this paper; 
this will be discussed in detail in the next section. It is worth pointing out that the incorporation of COI feedback into 
the hybrid- ARQ scheme is being implemented at the physical layer. After combining at the destination, the packet is 
then demodulated using either soft or hard decoding methods and then passed to the FEC decoder which then outputs 
a final estimate of the original information packet, m[fc]. 

It is important to note that a feedback channel is present between the destination and the source. In fact, for any 
ARQ protocol, a feedback channel is necessary so the destination can send back an ACK7NACK signal. In our setup, 
we assume that: 

• The destination does not only send back ACKVNACK information but also feeds back CSI which could be 
outdated, current, or quantized. 

• The destination can feed back the COI for the packet to the source where COI feedback is simply the destination 
feeding back exactly (or a subset of) what it has received. This is discussed further in Section V. 

Explicitly, the causal COI at the source is equivalent to the source having access to the past values of y[k]. However, 
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since noisy COI feedback is also investigated, we introduce a feedback noise process n[k] (see Fig. [TJ so that the source 
now only has access to past values of y[k] + n[k}. Note that the source might have access to all or only some of the 
entries in y [k] + n[fc] based on how much COI is being fed back. Furthermore, the source can subtract out what it sent 
due to the availability of h[k]. Therefore, this is analogous to having access to past values of z[k] + n[k]. The feedback 
noise, n[fc], is assumed to be complex AWGN such that n[k] ~ CAf(0, cr 2 I) and also independent of the forward noise 
process, z[k]. Note that setting a 2 = yields perfect COI feedback as a special case. 

III. Linear Feedback Combining 

We now narrow our focus to the packet encoding/combining steps of the hybrid-ARQ system; the full system 
including the FEC will be considered in Section V. Specifically, we consider employing COI feedback to better refine 
the destination's packet estimate 0[k] after each retransmission. Improving the quality of the estimate will lead to fewer 
decoding errors and higher throughput. To begin, we look at the most straightforward setup of SISO, where the source 
and destination each have one antenna and outdated CSI along with causal COI, whether it be noisy or perfect. The 
scheme will be extended for use with multiple antenna scenarios and the effects of varying the CSI feedback will be 
discussed in Section IV. 

A. Overview 

To construct the new transmission strategy, we aim to develop a linear coding scheme with the objective of 
maximizing the post-processed signal-to-noise ratio (SNR) after N retransmissions. Initially, we focus on sending 
only one symbol or, in other words, assume that L = 1 where the information packet is now a scalar, 9 G C. Note 
that this assumption is only done to reduce the amount of notation — the proposed scheme can be readily extended to 
arbitrary packet lengths as it is assumed that it will be utilized, in general, for L ^S> 1. In the case that L = 1, however, 
the transmit and received vectors x.[k] and y[fc] also reduce to scalars x[k] and y[k]. 

It is helpful at this stage to introduce a mathematical framework for linear feedback coding. It can be seen that if 
L = 1, then, gathering all packet transmissions together, (f2) can be rewritten as 

y = Dx + z, (3) 

where y = [y[l], y[2], . . . , y[iV]] T is a column vector (likewise for x and z) and D = diag(/i[l], h[2], . . . , h[N]) is a 
matrix formed with the channel coefficients down the diagonal. Note that the notation D is chosen to give distinction 
between it and the commonly-used H for a MIMO channel matrix which is used later in the paper. With this setup, 
we can write the transmit vector x as 

x = g 6» + F(z + n), (4) 

where g g C JVxl is the vector used to encode the symbol to be sent, 6, and F e i^nxn - ls a strictly lower triangular 
matrix used to encode the side-information {z + n}. The form of F is constrained to be strictly lower triangular to 
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enforce causality. Note that is the transmit structure of linear feedback coding — the transmitted value is a linear 
function of the side-information and of the information message. Furthermore, the encoding process is encapsulated 
by the matrix F and the vector g. 

Now, we shift focus to the destination. After receiving a packet, the destination forms an estimate of the packet by 

0[fc] = q*;y(fc): (5) 

where 8[k] is the destination's estimate of the symbol 9 after k retransmissions, qfe 6 C Nxl is the combining vector 
used for the k th retransmission, and the notation yvw refers to the first k entries of y. It can be seen that the packet 
estimation process can be completely described by the vectors q^. Thus, the entire linear feedback code can be 
represented by the tuple (g, F, q&) ll33l . Thus, with this framework, the post-processed SNR for the system after 
N retransmissions can be defined as 

o NR = lq* p gl 2 p rfi x 

||q*(I + DF)|| 2 + a 2 ||q*DF|| 2 ' W 
where q = qAr; the subscript is dropped for convenience. We can now mathematically define the overall objective of 
this section: to construct a linear feedback code that maximizes the post-processed SNR given the channel coefficients 
h[k]; this is equivalent to finding 

, v v |q*Dg|V n , 

(g, F, q) opt = (7) 

while satisfying the average power constraint ([T) and causality of side-information. 

In the special case of a SISO system with perfect COI (a 2 — 0) and causal CSI, it can be shown that the solution to 
(|7]i has the structure of the scheme given in ll35l (the specific derivation is omitted due to space concerns). Interestingly, 
the scheme in that work was derived using a control-theoretic approach instead of using post-processed SNR as an 
objective function ||36ll . Most importantly, the optimality of the scheme in the SNR sense motivates our construction. 
In particular, we develop a generalization of the feedback scheme presented in [35] as this scheme was not only shown 
to achieve capacity but also achieve a doubly exponential decay in probability of error. In the proposed generalized 
scheme, we extend the original scheme for use with: 

• single and multiple antennas (i.e., MISO and MIMO wireless systems), 

• perfect and noisy COI feedback, 

• outdated and current CSI at the source. 

The details of the proposed scheme are given in the following sections. 

B. Encoding 

The fundamental idea of the transmission scheme is to transmit the scaled estimation error from the previous 
transmission for each successive retransmission so that the destination can attempt to correct its current estimate ||3~7l . 
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|. The scaling operation is performed so the transmitted signal meets the average power constraint ([T}. To further 
illustrate this concept and help motivate our construction, we now briefly present a heuristic overview of the scheme 
in ll35l . In this case, the transmitted signal, x[k], is given as 

x[k + l] = 8[k]e[k], (8) 

where e[k] = 9 — 6[k] is the error in the destination's estimate of the message after the k th packet reception and 
S[k] is the scaling factor chosen to appease the power constraint. After receiving y[k], the destination then forms an 
estimate of the error, e[k]. This is then subtracted from the current estimate. As will be shown, our proposed scheme 
is motivated by this error-scaling technique. 

We now define how the source encodes the message. As the encoding process for a perfect COI feedback and 
the encoding process for a noisy COI feedback are very similar, we now introduce the encoding process for both 
perfect and noisy COI feedback in a single framework. The encoding operation of the proposed scheme can be written 
compactly in the definitions of F and g; they are constructed as: 

. The (i, j) th entry of F, foj, is 

-*fyp<f)[i- l]h*[j], i>j, 
0, i < j, 



fi,j 



• The i entry of g, gi, is 



1], 



where 



tf>[k\ 



J]> (7 , CT 2)[i], fc>0 

i=l 

1, k = Q, 



P h ,^)[k} = (1 + (l + a 2 ) 1P \h[k] | 2 ) 1/2 , (9) 

and 7 € [0, 1] is a constant. Note that the scaling factor S[k] in ([H} is now given its analog by the term <fi[k] which 
ensures the proposed scheme meets the power constraint ([T}. 

The scheme presented here in the form of g and F is a direct generalization of the error-scaling scheme in d8|) as the 
original scheme for perfect COI feedback can be obtained as a special case of these definitions by letting 7=1 and 
cr 2 = 0. The main mechanism introduced into the proposed scheme is a power allocation variable, 7, to help combat 
the effect of the feedback noise, n[k] J33j. Specifically, 7 is a degree of freedom introduced to allocate power between 
the encoding of feedback side-information and the information to be sent. It is only of use when the feedback channel 
is noisy; if feedback noise is not present, it should be set to 7 — 1 and disregarded. In brief, as 7 — > 0, this scheme 
simply repeats the packet on every retransmission (i.e., the scheme becomes equivalent to Chase combining). As 7 
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grows, the scheme uses most of the feedback power to mitigate the noise in the destination estimate. This quantity is 
discussed in detail later in the paper. 

Now, with g and F defined, the encoding process is completely described, and we can now move on to verifying that 
it meets the average transmit power constraint ((T). As will be shown, it is much easier to derive the average transmit 
power of the proposed scheme if it is rewritten in a recursive manner; thus, its recursive form is now presented. 
Assuming that the symbol is scaled such that £7[| 2 ] = p, the first packet transmission is set to the symbol itself with 
x[l] = 8. The subsequent transmissions can be written as 



x[k + 1] = /8( 7l0 a)[fc] {x[k} - y/jph*[k](z[k] + n[k])) , 1< k < N. 
With the recursive formulation given, we can now present the following lemma. 



(10) 



Lemma 1. The proposed scheme in f liOD meets the average transmit power constraint given in (fJJ for both noisy and 
perfect COI feedback. 

Proof: The proof is based on a simple inductive argument. Since the symbol has been scaled to have a second 
moment of p, the average power of the first transmission is | 2 ] = p. Assume that ^[^[fc]! 2 ] = p for some 

k. Using ( [Tol l, we can write the average transmit power for the (k + l) th retransmission of packet 9 conditioned on 
channel realization h[k] as 



E 



|i# + l]| 2 h[k 



= E 



|j8 (7 ,oa)[fc] (x[k] - ^fph*[k](z[k]+n[k}))\ h[k 
1 



(a) 



l + (l + (T 2 )7p|/l[fc]| 2 

1 

l + (l + (T 2 )7p|/l[fc]| 2 

1 

1 

p, 



E 



E 



x[k] - ^iph*[k}{z[k] +n[k])\ z h[k] 
:[k]\ 2 \h[k]\+(l + a 2 hp 2 \h[k}\ 2 ) 



{l + o*) 1P \h[k]f 



{p+{l + (j 2 hp 2 \h[kf 



where the equality in (a) follows from E \x[k]\ 2 h[k] — E [\x[k] | 2 ] = p. Therefore by the principle of mathematical 
induction, the equality holds for any arbitrary k. ■ 
Now that the encoding operation has been described and verified to meet the average transmit power constraint, it 
is possible to move on to the decoding stage. 

C. Decoding 

In this section, we discuss the decoding process in the proposed scheme. It is worthwhile to point out that we only 
perform soft signal-level decoding — the output of the destination is an estimate that is not necessarily mapped to an 
output alphabet. Unlike the encoding operation, decoding at the destination significantly differs depending on whether 
perfect or noisy COI is available at the source. 
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1) Perfect COI Decoding (a 1 = 0): First, we look into defining q for perfect COI. In the special case when 
the feedback channel is perfect, this scheme assumes the structure of the feedback scheme in [35|; we reproduce it in 
this section for completeness. In this case, the combining vector q has a concise closed form. In particular using the 
definition in (O, the i th component of q, q, L can be given as 

gi = #-lL8f li0) (ID 

Because of this definition, q/j for perfect COI can be defined as <\k — [qi, ■ ■ ■ , <Zfc] T - Note that since the COI at the 
source is assumed to be perfect, a 1 — and 7 = 1. Now that q has been defined, the entire scheme for perfect 
COI feedback has been described, and the structure of q can be used to formulate the decoding process in a recursive 
fashion. Thus, at this point, we introduce the following lemma: 

Lemma 2. The coding scheme for perfect COI feedback can be alternatively represented as 

x[k + l] = I3 {m [k}(x[k}- ph*[k}z[k}) (12) 

fc 

8[k] = (l-\4>[k]\ 2 ) 6 + p\(p[k]\ 2 ^ ^[m - l]h*[m]z[m\. (13) 

m— 1 

The proof has been relegated to the Appendix. Note that Lemma 2 suggests that the estimator of the proposed scheme 
is a biased one. However, we can easily make the final estimated output unbiased by performing the appropriate 
scaling. We can define the unbiased estimator of packet 8 as 

fr[k] = (i-i^MIT 1 ^] 

k 

= 6 + p(l - |0[fc]| 2 ) _1 \(j)[k}\ 2 -TV" l}h*[m]z[m}. (14) 

m— 1 

2) Noisy COI Decoding (a 2 > 0): The source is now assumed to have corrupted COI from the destination. 
Note that the two main differences between perfect COI decoding and noisy COI decoding are: 

• The power allocation variable, 7, is now a degree of freedom. This allows the source to allocate more or less 
power to the message signal to adapt to conditions of the feedback channel. 

• The destination can no longer be derived in a simple form as in the noiseless feedback case. It is derived from 
the form of the optimal linear estimator of the symbol, 8. 

It can be shown that, if a 1 > 0, the optimal q that maximizes post-processed SNR with the setup in (O and is 
given by 

C-^Dg .... 
g*D*C x Dg 

where C = (DF + I)(DF + 1)* + er 2 DFF*D* is the effective noise covariance matrix seen at the destination. Note 
that this definition of q assumes that all TV retransmissions are used as q will be a length N vector. To obtain q^ where 
1 < k < N, one can simply truncate the vectors (and matrices) in ([15} to simply the first k entries (rows and columns). 
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With this setup, the post-processed SNR, given the channel coefficients h[k], can be written as 

SNR = p (g*D*C _1 Dg) . (16) 

It is difficult to derive a simple expression for (fTol : we instead formulate bounds on the post-processed SNR. This is 
done in the following lemma for the case of N = 2 in the low and high p regimes. 

Lemma 3. Given the linear feedback code described above with blocklength N — 2, at small p (i.e., p <C 1), the 
average post-processed SNR can be bounded by 

E[SNR N=2 ]<2p(l + ^p + 1P 2 ), (17) 

and 

E[SNR N=2 ] > 2p fl + y/jp - —tt—Ip] ■ (18) 
Furthermore, at large p (i.e., p> 1), the average post-processed SNR expression behaves as: 

E[SNR N=2 ] — ► p ( 1 + -L V (19) 



Proof: In the case of TV = 2, the post-processed SNR using dTo"} can be calculated to be 

Using (O which states that /3( 7 , CT 2 ) [k] < 1 for any k, it is clear that, 

SNR N=2 < p (\h[l}\ 2 + \h[2]\ 2 (l + ^jp\h[l}ff) . 

Now taking expectation on both sides and using the independence of fading blocks in time and i?[|/i[l]| 2 ] =£'[|ft,[2]| 2 ] = 
1, E^l^l]! 4 ] = 2, we immediately get 

E[SNR N=2 ] < 2p (1 + V7P + IP 2 ) • (21) 
Using the inequality (1 + > (1 — £) valid for any real £ in d20i >. 

> p(|^[l]| 2 +/3^^ 2) [l]|^[2]| 2 (l + ^|/ l [l]| 2 ) 2 (l-cx^p^^^Cll^MI 2 ^^]! 2 )) . (22) 
Taking the conditional expectation with respect to h[2] in (1221 . we get 



E 



SNR N=2 



h[2]\ > p [\h[l}\ 2 + /3 2 7 ^)[1] (1 + /7HM1]| 2 ) (1 ~ 2a 2 7P 2 /3 2 7i(T2) [l]|Ml]| 2 ) 

= P (|Mi]| 2 + /3 2 7 , ff2) [i] (i + v^pIMi]! 2 ) 2 + o(p 2 )) . 

By the definition of /3 2 ^ a ^ [1] in (O and the inequality (1 + £) _1 > (1 — £), we have 

h[2]\ > p(|Ml]| 2 +(l-(l + ^ 2 )7p|^[l]| 2 )(l + /7p|^[l]| 2 ) 2 + 0(p 2 )) 
= p {\h[l}\ 2 + 1 + 2/7p|Ml]| 2 - (1 + ^hp\h[l}\ 2 + 0(p 2 )) . 



E 



SNR N=2 
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Now taking expectation with respect to the channel realization h[l], we immediately get 

E[SNR N=2 ] >2p(l + V7P (l - ^Y~V^j + 0(p 2 ) 
Therefore in the small p regime we have, 

E[SNR N=2 ] ^ 2p fl +^pfl- ^-^Vl 



Hence for the proposed linear scheme to have better performance than MRC, we require that 7 < , __ rr _. 
In the case of large p, the expression in d20l i by approximating the second term can be written as 

/3f 7(T2) [i]|M2]| 2 (i + V7HMi]l 2 ) 2 \ 



SNR N=2 — > p[\h[l]\' 



P [\h[i 



P [\h[i 



2 



2 



a 2 7P 2 /3 ( 2 7 , CT2) [l]|Ml]| 2 IM2]| 2 

(i + V7PlMi]l 2 ) 2 ^ 

a 2 jp 2 \h[l}\ 2 
1P 2 \h[l]\ i 



p-s-oo V <7 2 7/5 2 |/l[l]| 

= p(\h[l]\ 2 + ^\h[l]\- 



Now, taking expectation we get, 

E[SNR N=2 ] — > p (l + \ 
Note that when a 2 ss 0, the above scheme yields significant benefits. ■ 

D. Power Allocation 

In this subsection, we investigate the power allocation parameter 7 seen in the scheme for noisy COI feedback. 
As stated before, it can be roughly thought of as a measure of the amount of feedback side-information being used 
in the retransmission. Optimally choosing the value of 7 to maximize the post-processed SNR in ([Tol l is clearly a 
non-causal problem. Therefore instead we define ^ fadmg ) to be the one that maximizes the post-processed SNR over 
the ensemble average of all channel realizations, i.e., 

7 f ding) = max E [p(g*D*C- 1 Dg)] . (23) 

7G[0,1] 

The difficulty of analytically calculating the above quantity stems from the post-processed SNR having non-linear 
dependencies on the fading coefficients h[l], . . . , h[N]. However, it turns out that the optimal 7 in the i.i.d. Rayleigh 
block fading case (7p fadms ' ) ) is very close to the optimal 7 in the AWGN case (7q AWGN ' ) ) as derived in [33). This 
is displayed in Fig. [2] In Fig. [2] we see that the peaks of both performance curves for block fading (averaged over 
15,000 trials) and AWGN noise are quite close together. This is quite beneficial as it is very easy to numerically find 
the value of 7 that maximizes the post-processed SNR in the AWGN case, whereas it proves to be much more difficult 
in the presence of block fading. Because of the proximity of ^ AWGN ) anc j 7^ fadmg \ we assume that the value of 7 
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Fig. 2: Average post-processed SNR versus choice of power allocation, 7, for AWGN and block fading. 

that maximizes the average post-processed SNR, 70 = ^ AWGN ) ~ ^,( fadln s) _ T/he va i ue G f y Q does, however, change 
with the blocklength N. Furthermore, as the number of transmissions is not necessarily known ahead of time, it is 
intuitive to not choose 7 as a function of blocklength. Alternatively, we can fix 7 based on a reasonable number of 
packet retransmissions — this is discussed in the following example. 

Example 1 

To illustrate the performance of the linear feedback scheme, we now provide some simulations. In this first plot 
(Fig.Bat. the post-processed SNR of the scheme is plotted in contrast to MRC. MRC is analogous to using our scheme 
but setting 7 = 0. In other words, the source simply repeats the packet at each retransmission. Then, retransmissions 
are combined using a linear receiver similar to the one in (115) . The simulations were run with an average transmit 
power of p = 3 and for both noiseless COI feedback and varying levels of noisy COI feedback. As can be seen, the 
linear feedback outperforms MRC with a gap that increases with decreasing feedback noise. 

As mentioned above, 70 changes with blocklength, N, and therefore should be chosen appropriately. However, in 
hybrid-ARQ, the required number of retransmissions is often not known ahead of time. Despite this fact, not having 
this knowledge a priori provides very little penalty to performance. If the number of retransmissions is not assumed to 
be predetermined, 7 can be approximately chosen using the feedback noise variance a 2 and the average transmit power 
p. The next figure, Fig. [3b] illustrates the effect of fixing 7. As is illustrated, fixing 7 with respect to blocklength yields 
little performance degradation as long as 7 is chosen appropriately. The average post-processed SNR for 7 = 0.01 
performs very close to the scheme when using 70 from Fig. [2] Note that Fig. [3b] has been plotted on a linear scale to 
help display the comparison. 
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(a) Average Post-processed SNR (b) Power Allocation 

Fig. 3: (a) Average post-processed SNR performance of linear feedback combining (LFC) versus maximal 
ratio combining, (b) Average post-processed SNR performance of linear feedback combining using 
different values of 7. 



IV. Multiple Antenna Scenarios 

In this section, we show how the feedback scheme for SISO systems can be implemented in both MISO and MIMO 
systems with current CSI at the source. In addition, an extension of the scheme is given for MIMO systems with 
perfect COI and only outdated CSI at the source. However, first we look at a MISO system with current quantized CSI 
along with perfect COI available at the source. 

A. MISO with Current, Quantized Channel State Information at the Source 

Consider a MISO discrete-time system (Fig.lUi with Alt transmit antennas and only one receive antenna, where the 
received packet, y[fc] G C lxL is given by 

y[k] = h T [fc]X[fc] +z[k], k = l,...N, (24) 

where h[k] € C M * xl is the channel gain vector, X[fc] € c M *>< L j s the transmitted packet matrix where the columns 
correspond to channel uses and the rows correspond to antennas, and z[k] € C lxL is additive noise during the k th re- 
transmission with distribution CJV(0, 1). Furthermore, the power constraint at the source is given as E[tr(X* [fc]X[fc] )] < 
Lp, and it is assumed that there is perfect CSI at the destination. However, the source no longer has access to perfect 
CSI. The destination only feeds back the beamforming vector to be used for current packet retransmission. The 
previous channel quality information (h T [k — l]w[k — 1]) along with the unquantized channel output is also fed 
back to the source. 
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Fig. 4: System level block diagram for MISO system with quantized CSI and COI feedback. The feedback 
information can be described as COI (y[k — 1]), channel quality information (h T [k — l]w[k — 1]), and 
the beamforming vector (w[/c]). 



The transmitted packet matrix X[fc] is now generated as an outer product by 

X[fc] = w[fc]x[fc], (25) 

wherew[fc] € C MtXl denotes the unit norm beamforming vector to be used during retransmission fc and x[fc] € C lxL 
is the effective SISO signal during retransmission number fc. The power constraint on X[fc] now is equivalent to 

£[tr(X*[fc]X[fc])] = tr(J5[x*[fc]w*[fc]w[fc]x[fc]]) 

= mm] 2 } 

< Lp. (26) 

At this point, it is again assumed that L = 1 for simplicity which reduces y[fc], x[fc], and z[k] to scalars y[k], x[k], and 
z [fc] . We now follow the standard model for limited feedback beamforming by constraining the design of beamforming 
vector w[fc] for packet transmission fc to a codebook F\k\ containing 2 B unit vectors [f|. We denote the codebook 

J"[fc] as 

T[k] = {f 1 [k],...,f aB [k]}, ||f J [fc]|| = l,l<j<2 B . (27) 

We can use any scheme available in literature to generate the unit beamforming vectors including random vector 
quantization (RVQ) [39], ||40| and Grassmannian line packing BTI . Il42l . This codebook is accessible to both the 
source and destination simultaneously. For RVQ, there must be a random seed that is made available to both the 
source and destination before the communication starts. 

The destination decides on the beamforming vector that the source uses during the k th retransmission by solving 
the following channel quality maximization problem 

w[fc] = argmax |h T [fc]f i [fc]| 2 . (28) 

fj[k]eF[k] 

Effectively, the destination chooses the unit vector w[fc] in the codebook J-[k] along which the channel vector h[fc] has 
the largest projection. The information about w[fc] is conveyed back to the source in just B bits. The limited feedback 
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capacity (Clf) for a given codebook design {.F[fe]}fcLx can be expressed by 



max log 2 (l + p|h r [fc]f#]| 2 ) 



(29) 



log 2 (l+p max |h T [fc]f,[fc]| 2 ) 



Clf — ^ 

Using the monotonicity of the logarithmic function, Clf can be simplified to 

Clf = E 

= £[log 2 (l + /9 |h T [fc]w[fc]| 2 )] . (30) 

As the number of feedback bits B approach infinity, Clf — > CmisOi where Cmiso = E [log 2 (l + p||h|| 2 )] . This 
is because limited CSI feedback becomes perfect CSI feedback for any codebook design with an infinite number 
of feedback bits. With the selection of beamforming vector w[fc] as described above and packet length L = 1, the 
received signal y[k] is given as 

y[k] = h T [k}w[k}x[k] + z[k], k = l,...N. 

Pre-multiplying the received signal y[k] by e^J ZhT [' £ ] w [ fe ] ? we obtain 

y[k] = \h T [k]w[k]\x[k] +z[k], k = l,...N, 

where y [k] = y[k]e- jAh W W W and z[k] is distributed as CAf(0, 1). If we let X[k] = \h T [k]w[k]\, we get the overall 
system in d24b as 

y[k] = \[k]x[k] +z[k], k=l,...N, (31) 

Finally, gathering all packet retransmissions together as in (O, we can rewrite ( l3TT > as 

y = Ax + z, (32) 

where A = diag(A[l], A [2], . . . , A [TV]). With this formulation, the MISO system is equivalent to the SISO system in 
([3); therefore, the SISO scheme can be implemented by replacing the role of h[k] with A[fe] (see Fig.|4]i. 

It can be proven in a similar way to Lemma 1 that the MISO scheme meets the average power constraint. Also, 
it can be shown that if the feedback channel is perfect, the MISO scheme achieves the capacity of the channel and 
obtains a doubly exponential decay in error probability. However, to avoid redundancy, this proof is only given for 
the MIMO case in the next section (Lemma 4). The effects of using different vector quantization techniques and the 
overall performance of the MISO scheme are now presented in an example. 

Example 2 

To illustrate the potential of our scheme, consider a MISO system communicating over an i.i.d. Rayleigh block 
fading channel with each entry of h[fe] distributed as CAf(0, 1). In this example, the COI feedback is assumed to be 
noiseless (i.e., a 1 = 0). Using a limited CSI feedback framework, Fig. |5]plots the packet probability of error curves 
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Fig. 5: Variation of probability of error with the number of retransmissions for i.i.d. Rayleigh fading at 
p = dB and M t = 2. The performance of RVQ, Grassmanian line packing, and perfect feedback are 
compared for -R/Cmiso = 0.5 and R/Cmiso = 0.9 with B = 2 and B = 3. 
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Fig. 6: System level block diagram when channel state is known perfectly at both the source and 
destination. 



against the number of retransmissions for two different normalized rates of 0.5 and 0.9 where normalized rate is the 
ratio of the rate of transmission to the channel capacity (-R/Cmiso)- The plots are for p = OdB with a two-antenna 
source averaged over 10 6 i.i.d. fading realizations. The doubly exponential decay of the curves are clearly visible for 
all the feedback schemes: perfect CSI feedback and quantized CSI feedback - RVQ and Grassmanian line packing. 
Even with quantized CSI feedback and moderate normalized rate of 0.5, only a few retransmissions are required to 
achieve a very low packet error rate of 1% for both RVQ and Grassmanian line packing. 

B. MIMO with Current Channel State Information at the Source 

Consider now a MIMO packet retransmission system (Pig. |6j with M t transmit antennas and M r receive antennas 
where the number of spatial channels available is M = min(M r , M t ). The received matrix, Y[fc] £ <C MrXL , is given 
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by 

Y[jfe]=H[*]X[Jfe] + Z[Ar], (33) 

where X[fc] e C M * xi is, as in MISO, the transmit packet matrix, H[fc] e £M,,xM t is the y ck Rayleigh fading 
channel matrix whose entries are i.i.d. zero-mean complex Gaussian random variables with unit variance, and Z[k] g 
C ,x is an additive noise matrix with i.i.d. zero-mean complex Gaussian entries with unit variance. Note that due 
to the availability of multiple spatial channels, the total packet length has increased to contain ML symbols with M 
symbols transmitted over each channel use. Again, for the sake of simplicity, we assume that L = 1 which reduces 
Y[fc], X[fe], and Z[k] to column vectors y[k], x[fc], and z[k]. When the current block fading matrix is known both at 
the source and destination, we can effectively diagonalize the channel. Let 

H[fc] =U[fr]A[fc]V*[fc], (34) 

be a compact singular value decomposition (SVD) of the channel matrix H[fc], where \J[k] G C MrXtf , A[fc] £ 
C MxM , and V[k] e C M ' xM , with 

A[fc] = diag(Ai[fc],...,AM[fe]), \i[k]>\ 2 [k]...>\ M [k]>0 1 (35) 

U* [k]V[k] = V* [k] V[k] = I. (36) 

We can design the source vector x[k] as 

x[k]=\[k]s[k], (37) 

where s[k] € C Mxl with V[fe] defined by d34"i > and (l36b . Also pre-multiplying the received vector by U* [k], we obtain 
the effective system described by ( T33l > as 

U*[fc]y[fc] = U*[fc]H[fc]V[fc]s[jfc]+U*[fc]z[fe] 

y[k] = A[k]s[k] + z[k], (38) 

where y[k] G c Mxl anc [ z[k] E C . The effective noise z[k] is distributed as CAf(0,T) due to the rotational 
invariance of complex i.i.d. Gaussian vectors. Due to the a priori knowledge of the channel at the source, spatial 
waterfilling can be performed across the M parallel spatial channels for each packet transmitted. The entries of the 
waterfilling matrix S [k] = diag (£i [k] , . . . , £m [k] ) are defined as 

e#]=max(o ! ^-^), l<i<M. (39) 

The value of the constant £o [k] is the water-filling level chosen to satisfy the power constraint 

M 



£&[A] = 1. (40) 
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Furthermore, the capacity Ctr of a MIMO channel with the fading matrix known both at the source and destination 
can be written as 

M 

C T R = 5>[log 2 (l + ^A 4 2 )] , (41) 
1=1 

where we have dropped the retransmission index k due to the i.i.d. nature of the block fading matrix. With current 
CSI at the source and destination, the overall channel capacity of the MIMO channel can be expressed as a sum of M 
parallel non-interfering SISO spatial channels each with capacity C,; where Ci — E [log 2 (l + p£iA 2 )] , 1 < i < M. 
With the aid of the waterfilling matrix defined in ((39), (l38l can now be written as 

y[k] = A[fc]E[fc] 1/2 x]fc] +z[k], 

where s[k] = H[fc] 1 / 2 x[fc]. Note that the spatial waterfilling (or power adaptation) does not make use of the COI fed 
back to the source at all. Letting A[k] = A[fc]H[fc] 1 / 2 , the overall system can be represented in matrix form as 

y[k] = A[fc]x[fc] + z[k]. (42) 

We next transmit M symbols over M parallel spatial channels by exploiting the COI and previous CSI available at the 
source using a maximum of N transmissions. In other words, with (l42t . we can implement M parallel instances of 
the COI feedback SISO scheme — one for each spatial channel. Similar to the MISO case, we replace the role of h[k] 
with Xi[k] for the i th spatial channel. 

It is quite possible that each of the source constellations &i[N] has a different number of constellation points; 
note that &i[N] denotes the source constellation used for the i fh spatial channel. The number of equally likely 
constellation points chosen for each channel depends on the spatial capacity Cj of the subchannel. Therefore, the 
number of constellations points must be less than 2 NCi . 

The overall schematic of the proposed scheme, shown in Fig. [6] clearly demonstrates the independent constellation 
mapping of each of the M symbols of packet 6 along with the separation of the channel output adaptation from 
current channel state adaptation. Furthermore, it can be shown that, if the feedback channel is perfect, any rate less 
than capacity can be achieved by the above scheme at doubly exponential rate. 

Lemma 4. // a 2 = 0, the proposed scheme achieves any rate R < Ctr. Viewing the rate R as a sum of M 
spatial channel rates, R = X)f=i Rii tne coding scheme can achieve any rate Ri < Cifor the i th spatial channel. 
Furthermore the probability of error ( P e ) for the packet decays doubly exponentially as the function of the number of 
transmissions N. In other words, for sufficiently large N, 

Pe<Pi cxp (-2(' r * + W) , 
where fix and /3a are positive constants, while P3 is a real constant for a given rate R. 

Proof: See Appendix. ■ 
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C. MIMO with Outdated Channel State Information at the Source 

In the case when there are multiple antennas at both the source and destination and the source has access to only 
outdated CSI, a direct extension of the SISO scheme for perfect COI (LemmaO can be made. Using the same system 
setup as in ( |33l l, if L = 1, we can write the feedback scheme recursively as 

x[fc + l] = (1 + pH*[/s]H[fc]r 1/2 (x[k]~ P U*[k]z[k}) (43) 

k 

6[k] = (l-#[fc]**[jfc])0 + p#[fc]**[jfe] ^^^[m-lD'H'HzH, (44) 

m— 1 

where 

#[*] = (1 + MpU*[l]H[l}y 1/2 ■ ■ ■ (1 + MpU*[k]U[k]y 1/2 . (45) 

Unfortunately, even if the feedback channel is perfect but only outdated CSI is available at the source, it is difficult 
to prove a result similar to Lemma 4. Although it can be shown for some positive rates a doubly exponential decay of 
probability of error is achievable, it has not been proven for all rates below capacity. We now broaden our focus back 
to the view of the whole hybrid- ARQ scheme in the next section. 

V. The Hybrid-ARQ Scheme and Variations 

Rather than focusing on the packet estimate, 0[k], we now consider the overall hybrid- ARQ scheme including the 
FEC. For the FEC, we assume the use of a systematic turbo code; although, any systematic block code can work. 
Because of this choice, we will perform only soft decoding at the output of the packet combining step. This means for 
each symbol in the estimated packet, 9i [k] , we will form a set of log-likelihood ratios as 

LLR L = | log 

where ipj for j = 1,2, ... , \ Q[N] \ are the points of the constellation Q[N] utilized for modulation. Upon calculating 
these sets, they are then passed to the turbo decoder for decoding. The specific turbo code implemented is a rate 
Linfo/ Leaded = 1/3 turbo code defined in the UMTS standard; more details are given in the Simulations section. 

Now, we introduce different configurations of the overall scheme that might help adapt to different circumstances 
(e.g., feedback link rate, transmit/receive duration, etc.). To do so, we look at varying the amount of COI feedback 
sent to the source; this is also done to illustrate the trade-off between performance (e.g., throughput, FER, etc.) and the 
amount of information fed back. Note that the case of CSI-only feedback has already been explored in the literature; 
see for example ||2j, @, 1271 -1291. Therefore the emphasis here is in varying the amount of COI feedback. 

The most straightforward of the possible COI feedback configurations is one where the destination simply feeds 
back everything it receives without discrimination. This utilizes a noiseless/noisy version of the full received packet 
for feedback information; hence, we will refer to this method as full packet feedback (FPF). Alternatively, one can 
alter FPF by implementing a well-known concept in hybrid- ARQ with feedback (31); instead of feeding back all the 



p(0i = il>j\0i[k]) 



1,2,.. 



\Q[N]\ 



(46) 
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symbols of the received packet, we can instead feed back only the T "least reliable" symbols with their indices. The 
measure of "reliability" can be based off metrics such as the log-likelihood ratio (LLR) or the logarithm of the a 
posteriori probabilities (log-APP) [32 ] . Since only some of the symbols in the packet are fed back, we will refer to this 
scheme as partial packet feedback (PPF). 

A. Full Packet Feedback (FPF) 

In FPF, we propose a hybrid- ARQ scheme where the source is assumed to have access to a noiseless/noisy version 
of the last received packet. The performance of the perfect COI feedback scenario is used to demonstrate the maximum 
possible gains that can be achieved with the addition of channel output feedback. To help explicitly show the feedback 
information available, we now introduce r[k] as the channel output feedback side-information available at the source 
at packet transmission k. In FPF, 

r[k] =y[k- 1] + n[fc- 1], (47) 

where, in this case, r[k] G C lxL . 

As mentioned before, the first transmission of packet is assumed to be a codeword of a FEC code. If a NACK 
is received at the source, each subsequent packet is encoded symbol-wise by the linear feedback code described in 
Section III. This is used to refine the destination's estimate of each symbol in the original packet. To display the 
performance of the scheme, we look at comparing the normalized throughput of this scheme with the turbo-coded 
hybrid- ARQ used in 1 34 1 . This standard uses a rate-compatible punctured turbo code to encode the packet. Specifically, 
it uses a rate 1/3 UMTS turbo code ||431 and then punctures it for use in hybrid-ARQ. If sending one packet and 
M spatial channels are available for the MIMO setting, the assignment of M symbols for M spatial channels is 
done arbitrarily. Note that it is plausible that using dynamic adaptive modulation for each of the spatial channels or 
coordinating multiple retransmissions [2| might result in improvement in throughput. However, we do not consider 
this here, but we point out that in most of the cases our proposed scheme can be combined with the innovations on 
using CSI more efficiently. 

B. Partial Packet Feedback (PPF) 

For sake of practicality, it is desirable to minimize the amount of COI feedback information needed to be sent back 
to the source. As a step towards this, we now look at the effects of limiting the size of the COI feedback packet. 
We try to utilize the limited feedback channel in the most useful way by feeding back not the complete packet but 
only relatively few of the symbols in the received packet. As mentioned above, in the partial packet mode, the choice 
of COI feedback information is based on the relative reliability of soft decoded bits. This addition to the scheme is 
motivated by the technique used in [32 1 where it was shown that focusing on the least reliable information bits can 
greatly improve the performance of turbo-coded hybrid-ARQ. The selection process to construct the feedback packet, 
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r[k], is performed at the destination using the following method. The received packet y[k] (or Y[fc] as in d33l) ) is 
combined with the k — 1 previous received packets using MRC in the case of Chase combining or as described above 
if linear COI feedback coding is employed. After combining, the destination now has an estimate of the desired packet, 
0[k}. This packet estimate is now passed on to the turbo decoder, and its corresponding output is a set of LLRs for 
each original information bit. For notation, we refer to the LLR produced by the turbo decoder for the i th information 
bit, rrii, as £i, which can be mathematically written as 

r P K = o| y [i])" 



log 



1 < i < L inio . (48) 



jj(m. t = l|y[l])_ 

The least reliable bits are chosen as the T bits whose LLR values have the smallest magnitude, i.e., \£i\ < £ t h, 
where £ t h is chosen appropriately. The magnitude of ii close to zero indicates that the bit is almost equally likely to 
be either a 1 or a 0. Then, the set of symbols whose realizations are to be fed back is 

/ S ym - {0 k ■ rm € e k and < l th , 1 < i < L info , 1 < k < L} , (49) 

meaning the symbol is chosen to be fed back if it contains one of the least reliable information bits. Let |/ S ym | = T sym . 
With this technique, we can then write the feedback packet, rx[&], as 

r T [k] = y T [k-l] + n T [fc-l], (50) 

where 

y T [k - 1] = { Vi [k - 1] ■. 6i e I sym } . (51) 

Since T sym channel outputs are being fed back, n*r[fc] is of length T sym . Note that the selection process is straight- 
forward in this case because we assume the use of a systematic turbo code. It is also important to note that the 
Tsym symbols are chosen only once (after the first transmission). This process can be done after each retransmission 
but would require more feedback resources. Finally, if it is assumed that the number of channel uses per packet 
retransmission are constant, one can fill in the remaining L — T sym channel uses in numerous ways. One particular 
way is what we will refer to as partial packet feedback with partial Chase combining (PPF-PC). In this mode, on 
the forward transmission, the new T sym symbols generated for the T sym least reliable symbols based on our linear 
coding scheme are sent in conjunction with the repetition of the remaining (L — T sym ) other symbols used for Chase 
combining. 

VI. Simulations 

In this section we present numerical simulations to demonstrate the improvements possible with inclusion of our 
proposed linear COI feedback coding in hybrid-ARQ schemes. We assume that the channel is i.i.d. Rayleigh block 
fading. We limit the number of retransmissions to a maximum of four (i.e., N = 4). All the throughput calculations 
are done by averaging over 10 3 new packet transmissions. 
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The metric defined for calculating normalized throughput is given as: 

1 

where 1 < B < N is the number of transmissions needed for successful decoding of a packet 0. This can be 
equivalently thought of as a packet success rate or the inverse of the average number of packets needed for successful 
transmission. If the number of retransmissions reaches the maximum number before successful decoding, the through- 
put contribution is zero. Note that this metric is meaningful only when comparing constant-length packet schemes. 
Also the above throughput definition implies that as p — >• oo, r —> 1 for all the protocols; including Chase and our 
proposed scheme. 

Fig- EH compares the performance of FPF scheme with perfect COI feedback against Chase combining for QPSK, 
16-QAM, and 64-QAM constellations over a 2 x 2 MIMO channel. The FEC code used for simulations is a 1/3 UMTS 
turbo code with eight decoding iterations. It is seen that most of the gains from our proposed scheme are realized in 
the low SNR regime. The FPF for QPSK displays gains of around 1 dB over Chase combining, and in 16-QAM, it 
gives an improvement of about 2 dB over Chase combining. Furthermore, the gain increases to 3 dB when the denser 
constellation of 64-QAM is chosen. It should be noted that these gains have been realized directly at the packet level 
and not at the bit level. This shows that with four retransmissions the power required at the source can be halved with 
the inclusion of the proposed linear coding scheme. 

In Fig. [7b] we plot the normalized throughput for PPF-PC for perfect COI feedback against traditional Chase 
combining scheme for 16-QAM and 64-QAM over a SISO channel. The amount of COI feedback symbols from the 
destination to the source is varied from 33% to 75% of the total feedforward packet size. For 64-QAM with a 1/3 
UMTS turbo code and Lj n fo = 2020, the length of the packet is L = 1010 symbols. Therefore the number of feedback 
symbols T for 64-QAM is varied from T = 337 to T = 757. Again we can see the improvements for 16-QAM 
and 64-QAM. Although the gains are smaller than the ones for full packet feedback, they are still significant. It is 
actually interesting to note that in 16-QAM most of the improvement in performance is reached with only 50% of 
COI feedback information. With 33% COI feedback, PPF scheme still shows an improvement of 1 dB over Chase for 
16-QAM and a substantial improvement of 2 dB for 64-QAM constellation. 

Fig. |7c] compares the normalized throughput for FPF with noisy COI feedback against Chase combining and the 
scheme in [34|. It is seen that even with a noise of a 1 — 0.25 on the channel output feedback channel, we see 
an improvement of about 0.5 dB for the linear feedback scheme over the scheme in [34| in the low SNR regime. 
Furthermore this gain is realized with only the addition of a very low complexity linear coder at source and destination. 

Finally, Fig.|7d]compares the normalized throughput for FPF with quantized COI feedback for 16-QAM constella- 
tion. At the destination, each of the components - inphase and quardature phase, are quantized using certain number 
of bits and then mapped to the QAM constellation. It is seen that even with 1-bit of quantization per phase leads to 
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Fig. 7: (a) Plot of variation of throughput for FPF and Chase combining with channel SNR p for 2 x 2 
MIMO channel. The performance is compared for QPSK, 16-QAM and 64-QAM constellations with 
L in f M = 2020 bits, (b) Plot of variation of throughput for PPF and Chase combining with channel SNR 
p for a SISO channel. The amount of PPF feedback is varied from 33% to 75% of the total frame size. 
The performance is compared for 16-QAM and 64-QAM with -L in f = 2020 bits, (c) Plot of variation 
of throughput for FPF with noisy feedback against Chase combining and HSDPA for a SISO channel. 
The performance is compared for QPSK constellation with Li n fo = 3200 bits, (d) Plot of variation of 
throughput for FPF with quantized COI feedback against Chase combining for a SISO channel. The 
performance is compared for 16-QAM constellation with L in f = 2020 bits. 
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improvement in the performance over the conventional Chase combining. However, most of the gains are obtained for 
5-bits per phase of quantization. 

VII. Conclusions 

In this paper, we have investigated a new hybrid-ARQ scheme that utilizes COI feedback side-information from 
the destination. This is motivated by trying to close the performance gap between Chase combining and incremental 
redundancy using feedback in order to leverage the implementation savings of a Chase combining system ll44ll . In 
normal Chase combining, packets are combined using maximum ratio combining (MRC); however, the proposed 
scheme incorporates feedback by combining the packets using a linear feedback code for fading channels with noisy 
feedback. Note that this also includes a new encoding step. It was shown through Monte Carlo simulations that the 
post-processed SNR performance of the linear feedback scheme greatly outperforms that of regular MRC. In addition, 
since the code is built on linear operations, it adds little complexity to the overall packet encoder and decoder assuming 
feedback side information is present. The full hybrid-ARQ scheme was analyzed using two main modes of operation: 
full packet feedback (FPF) in which the source was assumed to have access to a noiseless/noisy version of the last 
received packet and partial packet feedback (PPF) in which only a subset of the received symbol are fed back to 
the source. Simulations show that the addition of feedback to hybrid-ARQ greatly increases the performance and 
outperforms incremental redundancy in most cases. 

Appendix 

A. Proof of Lemma |2] 

Proof: The encoding for perfect COI feedback using <j9j can be written for each x[k] as 

x[k + l] = -1 [fe]e[Ar], (52) 

where 

e[k] = 6-6[k]. (53) 
The operations at the decoder side can also be given by 

x[k + l] = (1 + p/i*[fc + l]/i[fc + l]) _1 ph*[k + l]y[k+l], (54) 

e[k] = 4>[k]x[k+ 1], (55) 

9[k + l] = 9[k] + e[k]. (56) 

For initialization purposes, it is assumed that 6[0] = 0. It can be seen from (|53l and (|56l l, that the error, e[k], for the 
symbol 9 satisfies the relation 

e[k + l] = e[k]-e[k]. (57) 
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Then, implementing (T52t and d57] i. we can rewrite x[k + 1] as 



According to 



x[k+l] = -1 [fc] (e[fe — 1] — e[fc — 1]) 

= (/T 1 [jfc] (<£[/c - l]x[k] - 4>[k - l]x[fc]) 

= {l + p\h[k)\ 2 ) X, \x[k]-x[k}) 

= (1 + pl^^l 2 )" 172 {x[k} - pfr*[fc]z[*c]) . 



0[fc] = (j)[k}x[k + 1] 

= tf,[k] (l + p\h[k]\ 2 y 1/2 (x[k]-ph*[k]z[k}) 
= cp[k] (1 + pl^]! 2 )- 172 x[k] - p<p[k] (1 + P \h[k]\ 2 y 1/2 h*[k]z[k] 

k 

= \4>[k]\ 2 Q - p\(j>[k]\ 2 <!>~ l \m - l]h*[m]z[m]. 

m— 1 

k 

6[k] = (l- \(t>[k)\ 2 ) 9 + p\<t>[k] | 2 ^ _1 [m- [m]z[m]. 

m— 1 

■ 

B. Proof of Lemma^ 

Proof: We present the proof for L = 1. The generalization of it immediately follows. For the i*' 1 spatial channel, 
we select the symbol 9i from a square QAM constellation consisting of Mj[iV] = 2 symbols. According to the 
recursive definition in Lemma 2, the i th spatial signal is given as 



Therefore, it follows that 



fc-i fc-i / fe-i . 

Xi[k] = ]J 6j ~pJ2 [U = ] ^H^N- (58) 



J=l \/l + pA 2 [/] m=l \(=m + p\ 2 [l 

Let 



Now ( 1581 ) can be rewritten as 

/ fe-i 



Based on (fl"4l i which describes the unbiased estimation algorithm at the destination, 



™ = 9 < + '7^£l¥V* (61) 
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Let 



i[N] 



Given channel realizations over blocklength N, {H[fc]}^ =1 , and a known 8 i7 the random variable ef[N] is just a 
complex Gaussian random variable with conditional mean 



E 



e«[A] {H[k]}% =1 A 



E 
0. 



Similarly for the variance of e"[A], we obtain 



V i-^[-^ fc =i^[*-i] 



l [N] 



N 

E 



The symbol 0, is drawn from a square QAM constellation 0, [AT] given by, 

©47V] = ^[A] {±1 ± lj, ±1 ± 3j, , ± (0WipvI - l) ± (V^M - l) •?} . 

where the scaling factor a, [A] satisfies the power constraint at the source 

E[\6 i \ 2 ] = l(M i [N]-l)am=p. 



(62) 



(63) 



(64) 



A correct decision about 0, is made by the destination if the error e"[A] falls within the square (D^A]) of length 
2yfi$N\. Let 

P e ({H[fc]}f =1 ,^) = P (e«[A] £ □ 4 [A]|{H[fc]}f =1 ,^) . 

Clearly, 



+ 



P e ({H[fc]}f =1 ,^) < p(|$He( e r[A])|>y^V]|{H[fc]}f =1 ,^) 
p(pm(e?[A])| > v^[A]|{H[fc]}f =1 ,^) , 

where $Ke(e"[A]) and 3m(e"[A]) denote the real and imaginary part of ef[N] respectively. Using the identical 
distribution of the real and imaginary components of the error e"[A], we get 

/ , — \ 

P e ({H[fc]}f =1 ,^) <4Q 



ai[N] 



Clearly, 



' v \J Var($He( e -[A]) {H[fc]}f =1 ,^) 
Var(fRe(eV[JV])|{H[A]}j^ 1 ,(9 i ) = ^Var (e?[A] |{H[fc]}f =1 , ^ 
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Therefore, 



P e ({H[k]}» =1 , 



< 4Q 



IN 



3 1 



AT Af[fc] 



4Q( V / ^V]), 



where 



Oi[JV] = 



3 1 



(Mi[JV]-l)^[iV]Efc =1 



N AJ[fe] 



(65) 



(66) 



We next show that with probability 1, a»[JV] increases at least exponentially with N. From the definition of <j>i[N] in 
d59l we have < </>j[iV] < 1, iV > 0. Also the definition implies that the sequence {&[JV]}wLo is a monotonically 
decreasing sequence for arbitrary channel matrices. Hence by Theorem 3.14 in fl31l . the sequence {4>i[N]}^ =0 
converges. Also, 



, N 

E[log 2 fa[N}} = --£>[log 2 
fe=i 

= 

2 



i + A ? W 



Using d68| i and the strong law of large numbers (SLLN), we know that for any given e > 0, 3N\ such that 



P 



1 VJV>JVi. 



In particular, 



By the almost sure convergence of {(f>i[N]}^ =a to zero, we can choose N 2 such that 



P I 



1 VN > N 2 



Using SLLN again, we obtain that for a given e > 0, 3N% such that 



P (jZ A * 2 M <ViN(l + e)j =1 V7V>7V 3 , 



where 77^ = E[\f [k]]. Substituting the bounds given by d69ll, (T70b and d7TT > into the expression of a.;[iV] in 
obtain that V7V > ma,x{N 1 ,N 2 , N 3 } with probability 1, 



ai[N] > 



P 2 NR >2- N ( 1 -^r ]l N(l + e) 

X 2 N((l-e)C i -R i ) 
~ P T] t N(l + e) 



(67) 
(68) 



(69) 



(70) 



(71) 



we 
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The positive value e also satisfies the inequality, 

pr]iN(l + e) < 2 eNC \ ViV > N 4 . 

Clearly it follows that ViV > iV max 

a t [N]>2 N((1 - 2t)C ^ Ri \ 

where 7V max = max{ , N 2 , N 3 , N 4 }. 

Thus, we have shown that with probability one, the input parameter of the Q-function increases exponentially. 
Furthermore it is very well known that Q-function decays exponentially and can be bounded by, 

Q{x) < \z~ x2/2 i Vx > 0. 
From the above two equations we immediately obtain, 

P e ({H[ fc ]}f =1 ,^)<2e-^ ((1 - 2 ^-^. 

Note that we can choose e arbitrarily. Picking e < min M I (l — I? 1 ) guarantees that the decay is doubly expo- 
nential. ■ 
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